Two-color multistep cascading and parametric soliton-induced waveguides 
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We introduce the concept of two-color multistep cascading for vectorial parametric wave mixing 
in optical media with quadratic (second-order or x^) nonlinear response. We demonstrate that 
the multistep cascading allows light-guiding-light effects with quadratic spatial solitons. With the 
help of the so-called 'almost exact' analytical solutions, we describe the properties of parametric 
waveguides created by two-wave quadratic solitons. 
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Recent progress in the study of cascading effects in 
optical materials with quadratic (second-order or x^) 
nonlinear response has offered new opportunities for all- 
optical processing, optical communications, and optical 
solitons |l],|2| . Most of the studies of cascading effects 
employ parametric wave mixing processes with a sin- 
gle phase- matching and, as a result, two-step cascad- 
ing. For example, the two-step cascading associated with 
type I second-harmonic generation (SHG) includes the 
generation of the second harmonic (to + uj = 2uj) fol- 
lowed by reconstruction of the fundamental wave through 
the down-conversion frequency mixing (DFM) process 
(2lu — lu = uj). These two processes are governed by one 
phase-matched interaction and they differ only in the di- 
rection of power conversion. 

The idea to explore more than one simultaneous 
nearly phase-matched process, or double-phase-matched 
(DPM) wave interaction, became attractive only re- 
cently ||fi}, for the purposes of all-optical transistors, 
enhanced nonlinearity-induced phase shifts, and polar- 
ization switching. In particular, it was shown [Q that 
multistep cascading can be achieved by two second-order 
nonlinear cascading processes, SHG and sum-frequency 
mixing (SFM), and these two processes can also support 
a novel class of multi-color parametric solitons ||. The 
physics involved into the multistep cascading can be un- 
derstood by analyzing a chain of parametric processes: 
SHG {uj + uj = 2uj) -► SFM {uj + 2uj = 3uj) -► DFM 
(3w - uj = 2u>) -> DFM {2u - uj = uj). The main disad- 
vantage of this kind of parametric processes for applica- 
tions is that it requires nonlinear media transparent up 
to the third harmonic frequency. 

Then, the important question is: Can we find paramet- 
ric processes which involve only two frequencies but allow 
to get all advantages of multistep cascading ? In this 
Rapid Communication, we answer positively this ques- 
tion introducing the concept of two-color multistep cas- 
cading. We demonstrate a number of unique features of 
multistep parametric wave mixing which do not exist for 
the conventional two-step cascading. In particular, using 
one of the processes of two-color multistep cascading, we 



show how to introduce and explore the concept of light 
guiding light for quadratic spatial solitons, that has been 
analyzed earlier for Kerr-like spatial solitary waves j^] 
but seemed impossible for parametric interactions. For 
the first time to our knowledge, we find 'almost exact' 
analytical solutions for two-wave quadratic solitons and 
investigate, analytically and numerically, the properties 
of parametric waveguides created by quadratic spatial 
solitons in x' 2 ' nonlinear media. 

To introduce more than one parametric process involv- 
ing only two frequencies, we consider vectorial interac- 
tion of waves with different polarization. We denote two 
orthogonal polarization components of the fundamental 
frequency (FF) wave {uj\ — uj) as A and B, and two or- 
thogonal polarizations of the second harmonic (SH) wave 
{u>2 = 2uj), as S and T. Then, a simple multistep cascad- 
ing process consists of the following steps. First, the FF 
wave A generates the SH wave S via type I SHG pro- 
cess. Then, by down-conversion SA-B, the orthogonal 
FF wave B is generated. At last, the initial FF wave 
A is reconstructed by the processes SB-A or AB-S, SA- 
A. Two principal second-order processes AA-S and AB-S 
correspond to two different components of the x^ sus ~ 
ceptibility tensor, thus introducing additional degrees of 
freedom into the parametric interaction. 

TABLE I. Possible multistep cascading processes 





Principal 


Equivalent 


(a) 


(AA-S, AB-S) 


(BB-S, AB-S); (AA-T, AB-T) 
(BB-T, AB-T) 


(b) 


(AA-S, AB-T) 


(BB-S, AB-T); (AA-T, AB-S) 
(BB-T, AB-S) 


(c) 


(AA-S, BB-S) 


(AA-T, BB-T) 


(d) 


(AA-S, AA-T) 


(BB-S, BB-T) 



Different types of multistep cascading processes are sum- 
marized in Table Q. The processes in the row (a) of Table || 
described above and the multistep cascading introduced 
in Ref. j| are qualitatively similar, but the latter involves 
a third-harmonic wave. 
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To demonstrate some of the unique properties of the 
multistep cascading, we discuss here how it can be em- 
ployed for light-guiding-light effects in quadratic media. 
For this purpose, we consider the principal DPM process 
(c) (see Table |) in the planar slab-waveguide geome- 
try. Using the slowly varying envelope approximation 
with the assumption of zero absorption of all interacting 
waves, we obtain 



Oz ox z 
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Alkl J~z + 15 + + 2X2i?V Afe2Z = 0, 



(1) 



where xi,2 = 2fci<7i i 2, the nonlinear coupling coefficients 
at are proportional to the elements of the second-order 
susceptibility tensor, and Aki and A&2 are the corre- 
sponding wave- vector mismatch parameters. 

To simplify the system (Q), we look for its station- 
ary solutions and introduce the normalized envelopes 
u, v, and w according to the following relations, A — 
7iw enp(i/3z — |Afciz), B = 72V exp(i/3z — ^AA^-z), 
and S = 73W exp(2i/3z), where 7-f 1 = 2\\x\ 1 j^ 1 = 
2xq(xiX2) 1 ^ 2 , and 7^ = Xi x o> an d the longitudinal 
and transverse coordinates are measured in the units of 
zq = (f3 — Afci/2) -1 and xq = (zo/2ki) 1 ^ 2 , respectively. 
Then, we obtain a system of normalized equations, 



u w 



0. 



.du d 2 u 

dz dx 2 

.dv d 2 v 
— r tt^7 — a i v + X v w = 0, 



dz 
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dx 2 
d 2 w 
dx 2 



(2) 
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v 2 )=0, 



where X = (X2/X1), «i = W~ Afc 2 /2)(/3- Afci/2)" 1 , and 
a = 4y3 (/? — Afei/2) _1 . Equations (||) are the fundamen- 
tal model for describing any type of multistep cascading 
processes of the type (c) (see Table |). 

First of all, we notice that for v = (or, similarly, 
u = 0), the dimensionless model (^) coincides with the 
corresponding model for the two-step cascading due to 
type I SHG discussed earlier Jl],|2| , and its stationary so- 
lutions are defined by the equations for real u and w, 



dx 2 
d 2 w 
dx 2 



— u + uw = 0, 



(3) 



-u 2 = 0, 



that possess a one-parameter family of two-wave local- 
ized solutions (uq , wq ) found earlier numerically for any 
q / 1, and also known analytically for a = 1, uq(x) = 
(3/\/2)sech 2 (x/2) = y/2w (x) (see Ref. §). 



Then, in the small-amplitude approximation, the equa- 
tion for real orthogonally polarized FF wave v can be 
treated as an eigenvalue problem for an effective waveg- 
uide created by the SH field wq(x), 



d v 

TJ + IXWo(x) - u\\v = 0. 



(4) 



Therefore, an additional parametric process allows to 
propagate a probe beam of one polarization in an effec- 
tive waveguide created by a two-wave spatial soliton in 
a quadratic medium with FF component of another po- 
larization. However, this type of waveguide is different 
from what has been studied for Kerr-like solitons because 
it is coupled parametrically to the guided modes and, as a 
result, the physical picture of the guided modes is valid, 
rigorously speaking, only in the case of stationary phase- 
matched beams. As a result, the stability of the corre- 
sponding waveguide and localized modes of the orthog- 
onal polarization it guides is a key issue. In particular, 
the waveguide itself (i.e. two-wave parametric soliton) 
becomes unstable for a < a cx ~ 0.2 0. 

In order to find the guided modes of the parametric 
waveguide created by a two-wave quadratic soliton, wc 
have to solve Eq. (0) where the solution wq(x) is known 
numerically only. These solutions have been also de- 
scribed by the variational method (§], but the different 
types of the variational ansatz used do not provide a very 
good approximation for the soliton profile at all a. For 
our eigenvalue problem the function wq(x) defines 
parameters of the guided modes and, in order to obtain 
accurate results, it should be calculated as close as pos- 
sible to the exact solutions found numerically. To resolve 
this difficulty, below we suggest a novel 'almost exact' 
solution that would allow to solve analytically many of 
the problems involving guadratic solitons, including the 
eigenvalue problem (Q). 

First, we notice that from the exact result at a = 1 
and the asymptotic result for large a, w w u 2 /(2a), it 
follows that the SH component wq(x) of Eqs. (|J) remains 
almost self-similar for a > 1. Thus, we look for the SH 
field in the form wq(x) — w m sech 2 (x/p), where w m and 
p are unknown yet parameters. The solution for Uq(x) 
should be consistent with this choice of the shape for 
SH, and it is defined by the first (linear for u) equation 
of the system @. Therefore, we can take u in the form 
of the lowest guided mode, Uo(x) = u m sech p (x/p), that 
corresponds to an effective waveguide wq(x). By match- 
ing the asymptotics of these trial functions with those 
defined directly from Eqs. (^) at small and large x, we 
obtain the following solution, 

uq(x) = u m sedi p (x/p), wq(x) — w m sech 2 (x / p) , (5) 



(w m - 1) ' 



4(w m - ir 

(2 - w m ) 
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where all parameters are functions of a only. It it easy to 
verify that, for a cr < a < oo, the SH amplitude varies in 
the region 1.3 < w rn < 2, so that all the terms in Eq. (||) 
remain positive. 

It is really amazing that the analytical solution (||), 
(^) provides an excellent approximation for the profiles 
of the two-wave parametric solitons found numerically. 
Figures |l](a,b) show a comparison between the maximum 
amplitudes of the FF and SH components and selected 
soliton profiles, respectively. As a matter of fact, the nu- 
merical and analytical results on these plots are not dis- 
tinguishable, and that is why we show them differently, 
by continuous curves and crosses. For a < 1, the SH 
profile changes, but in the region a > a CT the approx- 
imate analytical solution is still very close to the exact 
numerical one: a relative error is less than 1%, for the 
amplitudes, and it does not exceed 3%, for the power 
components. That is why we define the analytical solu- 
tion given by Eqs. (jjj), (||) as 'almost exact'. Details of 
the derivation, as well as the analysis of the case a < 1, 
will be presented elsewhere H. 
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FIG. 1. Comparison between the numerical (continu- 
ous curves) and 'almost exact' analytical (crosses) solutions 
for two-wave (FF - dotted, SH - solid) parametric solitons: 
(a) Maximum amplitudes; (b) Two-wave soliton profile at 
a = 4. 

Now, the eigenvalue problem (Q) can be readily solved 
analytically. The eigenmode cutoff values are defined 
by the parameter a\ that takes one of the discrete val- 
ues, a[ n) = (s- nf/p 2 , where s = -(1/2) + [(1/4) + 
WmXP 2 ] 1 / 2 ■ Number n stands for the mode order (n — 
0, 1, . . .), and the localized solutions are possible provided 
ri < s. The profiles of the guided modes can be found 
analytically in the form 

v n (x) = Vsech s ~ n (x/p)H(-n, 2s - n + 1, s - n + 1; (), 

where £ = h[l — tanh(x/p)], V is the mode amplitude, 
and H is the hypergeometric function. 

According to these results, a two- wave parametric soli- 
ton creates, in general, a multi-mode waveguide and 
larger number of the guided modes is observed for smaller 
a. Figures |](a,b) show the dependence of the mode cutoff 



values a ™ vs. a, at fixed Xi ancl vs - t ne parameter x, at 
fixed a, respectively. For the case x = F the dependence 
has a simple form: = [1 — n(w m — l)] 2 . 




FIG. 2. Cutoff eigenvalues ot^ of the guided modes shown 
as (a) functions of a at \ = 2, and (b) functions of \ a t a = 4. 
Dashed lines correspond to the intersection of the plots in the 
parameter space (a,\)- 

Because a two-wave soliton creates an induced waveg- 
uide parametrically coupled to the modes of the orthog- 
onal polarization it guides, the dynamics of the guided 
modes may differ drastically from that of conventional 
waveguides based on the Kerr- type nonlinearities. Fig- 
ures |](a-d) show two examples of the evolution of guided 
modes. In the first example [see Fig. ||(a-c)], a weak fun- 
damental mode is amplified via parametric interaction 
with a soliton waveguide, and the mode experiences a 
strong power exchange with the orthogonally polarized 
FF component through the SH field, but with only a 
weak deformation of the induced waveguide [see Fig. |^(a) 
- dotted curve] . This effect can be interpreted as a power 
exchange between two guided modes of orthogonal polar- 
izations in a waveguide created by the SH field. In the 
second example, the propagation is stable [see Fig. §(d)]. 

When all the fields in Eq. (|||) are not small, i.e. 
the small-amplitude approximation is no longer valid, 
the profiles of the three-component solitons should be 
found numerically. However, some of the lowest-order 
states can be calculated approximately using the ap- 
proach of the 'almost exact' solution (|^),(||) described 
above. Moreover, a number of the solutions and their 
families can be obtained in an explicit analytical form. 
For example, for a\ — 1/4, there exist two families of 
three- component solitary waves for any a > 1, that de- 
scribe soliton branches starting at the bifurcation points 
on = c^p at a = 1: (i) the soliton with a zero-order 
guided mode for x = 1/3: u(x) = (3/\/2~) sech 2 (x/2), 
v{x) = c 2 sech (x/2), w{x) = (3/2) sech 2 (x/2), and 
(ii) the soliton with a first-order guided mode for % = 1: 
u(x) = c\ sech 2 (x/2), v(x) = c 2 sech 2 (x/2) sinh (x/2), 
w(x) = (3/2) sech 2 (x/2), where C2 = y/3 (a — 1) and 
ci = \J (9/2) + c\. Some other soliton solutions ex- 
ist for a specific choice of the parameters, e.g. for 
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a = ol\ = 4/9 and \ — 1, we find u(x) = 
(4/3) sech 3 (x/3), v(x) = (4/3) sech 3 (x/3) sinh (x/3), 
and w(x) = (4/3) sech 2 (x/3). Stability of these three- 
wave solitons is a nontrivial issue; a rigorous analysis of 
all such multi-component states is beyond the scope of 
the present Rapid Communication and will be addressed 
elsewhere. 

At last, we would like to mention that in the limit 
of large a, when the coupling to the second harmonic 
is weak, we can use the cascading approximation w ~ 
(u 2 + v 2 ) /(2a). Then, the equations for two orthogonal 
polarizations of the FF wave reduce to a system of two 
coupled NLS equations, an asymmetric case of TE-TM 
vector spatial solitons well studied in the literature (see, 
e.g., Ref. |[o| and references therein). 




FIG. 3. (a) Change of the normalized power in FF (u, 
solid) and SH (w, dotted) components, which initially consti- 
tute a two-wave soliton, and in the guided mode (v, dashed) 
at x = 2, demonstrating amplification of a guided wave. Evo- 
lution of the guided wave and effective waveguide (SH) is 
presented in plots (b) and (c), respectively, (d) Stationary 
propagation of a stable fundamental mode (\ — 1). For all 
the plots a = 4, the initial amplitude is Vo = 0.1, and ai 
corresponds to the bifurcation point. 

For a practical realization of the DPM processes and 
the soliton light-guiding-light effects described above, we 
can suggest two general methods. The first method is 
based on the use of two commensurable periods of the 
quasi-phase-matched (QPM) periodic grating. Indeed, 
to achieve DPM, we can employ the first-order QPM for 
one parametric process, and the third-order QPM, for 
the other parametric process. Taking, as an example, 
the parameters for LiNbC>3 and AA-S (xx — z) and BB-S 
(zz — z) processes pdj , we find two points for DPM at 



about 0.89 /im and 1.25 /mi. This means that a single 
QPM grating can provide simultaneous phase-matching 
for two parametric processes. For such a configuration, 
we obtain \ fz 1.92 or, interchanging the polarization 
components, x ~ 0.52. The second method to achieve 
the conditions of DPM processes is based on the idea of 
quasi-periodic QPM grating. As has been recently shown 
experimentally fll2] ] and numerically |i~3| ], Fibonacci op- 
tical superlattices provide an effective way to achieve 
phase-matching at several incommensurable periods al- 
lowing multi-frequency harmonic generation in a single 
structure. 

In conclusion, we have introduced the concept of two- 
color multistep cascading and demonstrated a possibility 
of light-guiding-light effects with parametric waveguides 
created by two-wave spatial solitons in quadratic media. 
We believe our results open a new direction in research 
of cascading effects, and may bring new ideas into other 
fields of nonlinear physics where parametric wave inter- 
actions are important. 
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